Imbibition of Femtoliter-Scale DNA-Rich Aqueous Droplets into Porous Nylon Substrates by Molecular Printing.

This work presents the first reported imbibition mechanism of femtoliter (fL)-scale droplets produced by microchannel cantilever spotting (µCS) of DNA molecular inks into porous substrates (hydrophilic nylon). Differently from macroscopic or picoliter droplets, the downscaling to the fL-size leads to an imbibition process controlled by the subtle interplay of evaporation, spreading, viscosity, and capillarity, with gravitational forces being quasi-negligible. In particular, the minimization of droplet evaporation, surface tension, and viscosity allows for a reproducible droplet imbibition process. The dwell time on the nylon surface permits further tuning of the droplet lateral size, in accord with liquid ink diffusion mechanisms. The functionality of the printed DNA molecules is demonstrated at different imbibed oligonucleotide concentrations by hybridization with a fluorolabeled complementary sequence, resulting in a homogeneous coverage of DNA within the imbibed droplet. This study represents a first step toward the µCS-enabled fabrication of DNA-based biosensors and microarrays into porous substrates.

The p53 tetramer shows an induced-fit interaction of the C-terminal domain with the DNA-binding domain.

The Trp53 gene is the most frequently mutated gene in all human cancers. Its protein product p53 is a very powerful transcription factor that can activate different biochemical pathways and affect the regulation of metabolism, senescence, DNA damage response, cell cycle and cell death. The understanding of its function at the molecular level could be of pivotal relevance for therapy. Investigation of long-range intra- and interdomain communications in the p53 tetramer-DNA complex was performed by means of an atomistic model that included the tetramerization helices in the C-terminal domain, the DNA-binding domains and a consensus DNA-binding site of 18 base pairs. Nonsymmetric dynamics are illustrated in the four DNA-binding domains, with loop L1 switching from inward to outward conformations with respect to the DNA major groove. Direct intra- and intermonomeric long-range communications between the tetramerization and DNA-binding domains are noted. These long-distance conformational changes link the C terminus with the DNA-binding domain and provide a biophysical rationale for the reported functional regulation of the p53 C-terminal region. A fine characterization of the DNA deformation caused by p53 binding is obtained, with 'static' deformations always present and measured by the slide parameter in the central thymine-adenine base pairs; we also detect 'dynamic' deformations switched on and off by particular p53 tetrameric conformations and measured by the roll and twist parameters in the same base pairs. These different conformations can indeed modulate the electrostatic potential isosurfaces of the whole p53-DNA complex. These results provide a molecular/biophysical understanding of the evident role of the C terminus in post-translational modification that regulates the transcriptional function of p53. Furthermore, the unstructured C terminus is able to facilitate contacts between the core DNA-binding domains of the tetramer.

Oxidovanadium(IV) complexes with chrysin and silibinin: anticancer activity and mechanisms of action in a human colon adenocarcinoma model.

Vanadium compounds were studied during recent years to be considered as a representative of a new class of nonplatinum metal antitumor agents in combination to its low toxicity. On the other hand, flavonoids are a wide family of polyphenolic compounds synthesized by plants that display many interesting biological effects. Since coordination of ligands to metals can improve the pharmacological properties, we report herein, for the first time, a exhaustive study of the mechanisms of action of two oxidovanadium(IV) complexes with the flavonoids: silibinin Na2[VO(silibinin)22]·6H2O (VOsil) and chrysin [VO(chrysin)2EtOH]2(VOchrys) on human colon adenocarcinoma derived cell line HT-29. The complexes inhibited the cell viability of colon adenocarcinoma cells in a dose dependent manner with a greater potency than that the free ligands and free metal, demonstrating the benefit of complexation. The decrease of the ratio of the amount of reduced glutathione to the amount of oxidized glutathione were involved in the deleterious effects of both complexes. Besides, VOchrys caused cell cycle arrest in G2/M phase while VOsil activated caspase 3 and triggering the cells directly to apoptosis. Moreover, VOsil diminished the NF-kB activation via increasing the sensitivity of cells to apoptosis. On the other hand, VOsil inhibited the topoisomerase IB activity concluding that this is important target involved in the anticancer vanadium effects. As a whole, the results presented herein demonstrate that VOsil has a stronger deleterious action than VOchrys on HT-29 cells, whereby suggesting that Vosil is the potentially best candidate for future use in alternative anti-tumor treatments.

MD simulations of papillomavirus DNA-E2 protein complexes hints at a protein structural code for DNA deformation.

The structural dynamics of the DNA binding domains of the human papillomavirus strain 16 and the bovine papillomavirus strain 1, complexed with their DNA targets, has been investigated by modeling, molecular dynamics simulations, and nuclear magnetic resonance analysis. The simulations underline different dynamical features of the protein scaffolds and a different mechanical interaction of the two proteins with DNA. The two protein structures, although very similar, show differences in the relative mobility of secondary structure elements. Protein structural analyses, principal component analysis, and geometrical and energetic DNA analyses indicate that the two transcription factors utilize a different strategy in DNA recognition and deformation. Results show that the protein indirect DNA readout is not only addressable to the DNA molecule flexibility but it is finely tuned by the mechanical and dynamical properties of the protein scaffold involved in the interaction.

Molecular dynamics of the DNA-binding domain of the papillomavirus E2 transcriptional regulator uncover differential properties for DNA target accommodation.

Papillomaviruses are small DNA tumor viruses that infect mammalian hosts, with consequences from benign to cancerous lesions. The Early protein 2 is the master regulator for the virus life cycle, participating in gene transcription, DNA replication, and viral episome migration. All of these functions rely on primary target recognition by its dimeric DNA-binding domain. In this work, we performed molecular dynamics simulations in order to gain insights into the structural dynamics of the DNA-binding domains of two prototypic strains, human papillomavirus strain 16 and the bovine papillomavirus strain 1. The simulations underline different dynamic features in the two proteins. The human papillomavirus strain 16 domain displays a higher flexibility of the beta2-beta3 connecting loop in comparison with the bovine papillomavirus strain 1 domain, with a consequent effect on the DNA-binding helices, and thus on the modulation of DNA recognition. A compact beta-barrel is found in human papillomavirus strain 16, whereas the bovine papillomavirus strain 1 protein is characterized by a loose beta-barrel with a large number of cavities filled by water, which provides great flexibility. The rigidity of the human papillomavirus strain 16 beta-barrel prevents protein deformation, and, as a consequence, deformable spacers are the preferred targets in complex formation. In contrast, in bovine papillomavirus strain 1, a more deformable beta-barrel confers greater adaptability to the protein, allowing the binding of less flexible DNA regions. The flexibility data are confirmed by the experimental NMR S2 values, which are reproduced well by calculation. This feature may provide the protein with an ability to discriminate between spacer sequences. Clearly, the deformability required for the formation of the Early protein 2 C-terminal DNA-binding domain-DNA complexes of various types is based not only on the rigidity of the base sequences in the DNA spacers, but also on the intrinsic deformability properties of each domain.

Gene expression profile study in CFTR mutated bronchial cell lines.

Cystic fibrosis (CF) is caused by mutations in the gene encoding the cystic fibrosis conductance transmembrane regulator (CFTR). Symptoms are pancreatic insufficiency, chronic obstructive lung disease, liver disease, chronic sinusitis and infertility in male patients. The phenotypic variability may be explained only in part by the more than 1200 CFTR mutations, which are grouped into six different classes, according to their effect on the protein ranging from a severe (no synthesis or blocked processing) to mild mutation (altered conductance or reduced synthesis). However, it is now accepted that other genes (CF modifiers) influence the phenotypic spectrum of the disease. In order to identify CF modifier genes, we built a low-density home-made oligoarray containing 144 genes selected according to biochemical criteria and evaluated their expression in two CF bronchial epithelial cell lines (CuFi1 F508del/F508del; CuFi3 F508del/R553X). If we consider both cell lines, 38 genes (26.3%) show an altered expression pattern with a threshold > +/-1.5. Of these 38 genes, 12 are altered in CuFi1, and 26 in CuFi3. Some of these genes share the same expression pattern in both cell lines, while others have a different behaviour. These results were validated by a QRT-PCR assay (R2 CuFi1 = 0.81 and R2 CuFi3 = 0.91). These data could suggest that the presence of a class I allele (R553X) determines a more profound alteration of gene expression pattern than the presence of a class II allele (F508del). The identification of the genes altered by a specific CF mutation could lead to the development of a pharmacological approach specific for different CFTR genotypes.

Structural dynamics of the mitochondrial ADP/ATP carrier revealed by molecular dynamics simulation studies.

The mitochondrial adenosine diphosphate/adenosine triphosphate (ADP/ATP) carrier has been recently crystallized in complex with its specific inhibitor carboxyatractyloside (CATR). In the crystal structure, the six-transmembrane helix bundle that defines the nucleotide translocation pathway is closed on the matrix side due to sharp kinks in the odd-numbered helices. The closed conformation is further sealed by the loops protruding into the matrix that interact through an intricate network of charge-pairs. To gain insight into its structural dynamics we performed molecular dynamics (MD) simulation studies of the ADP/ATP carrier with and without its cocrystallized inhibitor. The two trajectories sampled a conformational space around two different configurations characterized by distinct salt-bridge networks with a significant shift from inter- to intrarepeat bonding on the matrix side in the absence of CATR. Analysis of the geometrical parameters defining the transmembrane helices showed that even-numbered helices can undergo a face rotation, whereas odd-numbered helices can undergo a change in the wobble angle with a conserved proline acting as molecular hinge. Our results provide new information on the dynamical properties of the ADP/ATP carrier and for the first time yield a detailed picture of a stable carrier conformation in absence of the inhibitor.

Experimental and simulative dissociation of dimeric Cu,Zn superoxide dismutase doubly mutated at the intersubunit surface.

The equilibrium properties of dimeric Photobacterium leiognathi Cu,Zn superoxide dismutase mutant bearing two negative charges in the amino acid clusters at the association interface has been studied, experimentally and computationally, and compared to those of the native enzyme. Pressure-dependent dissociation is observed for the mutant, as observed by the fluorescence shift of the unique tryptophan residue located at the intersubunit surface. The spectral shift occurs slowly, reaching a plateau after 15-20 min, and is fully reversible. Measurement of the degree of dissociation allows us to calculate the standard volume variation upon association and the dissociation constant at atmospheric pressure. On the other hand the native protein is undissociable at any pressure. In the simulative approach, the dissociation free energy has been calculated through the blue moon calculation method for the case of a multidimensional reaction coordinate, corrected for the rotational contribution within the semiclassical approximation for a free rigid-body rotor. The scheme permits to define a definite path for the rupture of the dimer and to calculate the effective force involved in the process. The calculated free energy difference is close to the experimental one, and the value obtained for the mutant is well below that obtained for the native protein, indicating that the theoretical reaction scheme is able to reproduce the experimental trend. Moreover, we find that, when the separation distance increases, the protein structure of the monomer is stable in line with the fast recovery of the original fluorescence properties after decompression, which excludes the presence of partly unfolded intermediates during the dimer-monomer transition.

Pre-discharge stress echocardiography and exercise ECG for risk stratification after uncomplicated acute myocardial infarction: results of the COSTAMI-II (cost of strategies after myocardial infarction) trial.

OBJECTIVE: To compare in a prospective, randomised, multicentre trial the relative merits of pre-discharge exercise ECG and early pharmacological stress echocardiography concerning risk stratification and costs of treating patients with uncomplicated acute myocardial infarction. DESIGN: 262 patients from six participating centres with a recent uncomplicated myocardial infarction were randomly assigned to early (day 3-5) pharmacological stress echocardiography (n = 132) or conventional pre-discharge (day 7-9) maximum symptom limited exercise ECG (n = 130). RESULTS: No complication occurred during either stress echocardiography or exercise ECG. At one year follow up there were 26 events (1 death, 5 non-fatal reinfarctions, 20 patients with unstable angina requiring hospitalisation) in patients randomly assigned to early stress echocardiography and 18 events (2 reinfarctions, 16 unstable angina requiring hospitalisation) in the group randomly assigned to exercise ECG (not significant). The negative predictive value was 92% for stress echocardiography and 88% for exercise ECG (not significant). Total costs of the two strategies were similar (not significant). CONCLUSION: Early pharmacological stress echocardiography and conventional pre-discharge symptom limited exercise ECG have similar clinical outcome and costs after uncomplicated infarction. Early pharmacological stress echocardiography should be considered a valid alternative even for patients with interpretable baseline ECG who can exercise.

Evidence of domain formation in cardiolipin-glycerophospholipid mixed monolayers. A thermodynamic and AFM study.

The interaction of the three main components of the mitochondrial membrane, namely cardiolipin, phosphatidylcholine, and phosphatidylethanolamine, has been studied investigating mixed cardiolipin-phosphatidylcholine and cardiolipin-phosphatidylethanolamine monolayers at different cardiolipin molar fractions. The thermodynamic behavior of the mixed monolayers was investigated by means of surface pressure and surface potential measurements, and atomic force microscopy was employed to characterize the morphology of the monolayers. Langmuir isotherms and surface potential curves show a regular behavior with a progressive transition toward the isotherm of the pure component. Positive deviations from ideality in the excess Gibbs energies of mixing suggest the presence of repulsive interactions in both systems. Analysis of partial molecular dipole moment indicates a discontinuity at a definite cardiolipin/phosphatidylethanolamine molar fraction, suggesting the formation of a stoichiometric complex; as a consequence, in mixed cardiolipin-phosphatidylethanolamine monolayers, a phase separation is observed at phosphatidylethanolamine excess. AFM measurements indicate the presence of two domains: one made by phosphatidylethanolamine and the other by a regular arrangement of phosphatidylethanolamine and cardiolipin at a fixed molecular ratio.

Role of the linker domain and the 203-214 N-terminal residues in the human topoisomerase I DNA complex dynamics.

The influence of the N-terminal residues 203-214 and the linker domain on motions in the human topoisomerase I-DNA complex has been investigated by comparing the molecular dynamics simulations of the system with (topo70) or without (topo58/6.3) these regions. Topo58/6.3 is found to fluctuate more than topo70, indicating that the presence of the N-terminal residues and the linker domain dampen the core and C-terminal fluctuations. The simulations also show that residues 203-207 and the linker domain participate in a network of correlated movements with key regions of the enzyme, involved in the human topoisomerase I catalytic cycle, providing a structural-dynamical explanation for the better DNA relaxation activity of topo70 when compared to topo58/6.3. The data have been examined in relation to a wealth of biochemical, site-directed mutagenesis and crystallographic data on human topoisomerase I. The simulations finally show the occurrence of a network of direct and water mediated hydrogen bonds in the proximity of the active site, and the presence of a water molecule in the appropriate position to accept a proton from the catalytic Tyr-723 residue, suggesting that water molecules have an important role in the stabilization and function of this enzyme.

Thermophilicity of wild type and mutant cold shock proteins by molecular dynamics simulation.

In the attempt to clarify possible mechanisms underlying thermal stability of proteins, we study through molecular dynamics thermophile Bc-Csp, mesophile Bs-CspB, and selected mutants. These proteins have been extensively characterized experimentally; researchers showed that differential thermostability among the wild type proteins is fundamentally linked to one or two mutated amino acids, and that the nature of the effect is electrostatic. They also inferred an atomistic mechanism related to removal of unfavorable interactions, rather than to the formation of salt bridges. Molecular dynamics allows us to confirm and support both hypotheses. Several other collective parameters have also been monitored in relation to thermophilicity, such as global and local rigidity, permanence and number of hydrogen bonds, or of salt links. None of these clearly correlates with the thermal stability of the presently studied proteins.

Prokaryotic Cu,Zn superoxide dismutases.

The Cu,ZnSODs (Cu,Zn superoxide dismutases) comprise a class of ubiquitous metalloenzymes that catalyse the dismutation of the superoxide radical anion into oxygen and hydrogen peroxide. The dismutation reaction involves two successive encounters of the superoxide anion with a catalytic copper centre hosted by the enzyme at the dead end of a narrow protein channel. Cu,ZnSOD is found in all eukaryotic species as a homodimeric enzyme of approx. 2x16 kDa. Cu,ZnSODs are also widely distributed in the prokaryotic phyla, and are found in the periplasmic space of bacterial pathogens, where they are likely to play a defensive role against extracellular oxidative stress, particularly with respect to phagocytic processes. The crystallographic structures of prokaryotic Cu,ZnSODs have shown that the core structural organization of the enzyme is based on a flattened antiparallel beta-barrel, composed of eight strands arranged in a Greek-key topological order, which usually forms a dimeric structure through the interaction of residues provided by beta-strands 4f and 5e with the 2,3 and 7,8 loops. The interface residues are arranged to form a closed ring, which surrounds a central cavity filled by water molecules that are arranged differently depending on the species. The interface variability of residues and solvent structure observed in bacterial enzymes permits fine tuning of the subunit association that results in the presence of monomeric and dimeric species, and confers a particular protein flexibility that gives rise to long-range communications between different region of the enzyme.

Prognostic and clinical correlates of angiographically diffuse non-obstructive coronary lesions.

OBJECTIVE: To make a prospective assessment of the clinical and prognostic correlates of angiographically diffuse non-obstructive coronary lesions. DESIGN: Angiographic vessel and extent scores were calculated in 228 clinically stable patients (mean (SD) age, 60 (11) years; 43 women, 185 men) undergoing prospective follow up for the composite end point of death and myocardial infarction. The effect on outcome of clinical variables (age, sex, previous myocardial infarction, diabetes mellitus, smoking habit, systemic hypertension, hypercholesterolaemia, ejection fraction) and angiographic variables (vessel and extent score) was evaluated by Cox's proportion hazard model. RESULTS: The vessel score was 3 in 34 patients (15%), 2 in 78 (34%), 1 in 87 (38%), and 0 in 29 (13%). Median extent score was 60 (range 6-110; first quartile 40, third quartile 70). Forty one events (nine deaths and 32 myocardial infarcts) occurred over a median follow up period of 30 months. Age and extent score were the only multivariate predictors of outcome, but the latter provided 28% additional prognostic information after adjustment for the most predictive variables (gain in chi2 = 7, p < 0.01). A vessel score of 3 was associated with worse survival, while no significant discrimination was possible among the other groups. However, assignment of patients to two groups according to an ROC curve derived cut off value for the extent score made it possible to obtain significant discrimination of survival even in cases with vessel scores of 0 to 2. Age and diabetes were clinical markers of a higher extent score. CONCLUSIONS: The angiographic extent score is a powerful marker of adverse outcome independent of severity and the number of flow limiting coronary lesions, and may reflect the link between clinical risk profile and diffusion of coronary atherosclerosis. Thus it should be of clinical value for targeting aggressive preventive measures.

Static and dynamic water molecules in Cu,Zn superoxide dismutase.

Understanding protein hydration is a crucial, and often underestimated issue, in unraveling protein function. Molecular dynamics (MD) computer simulation can provide a microscopic description of the water behavior. We have applied such a simulative approach to dimeric Photobacterium leiognathi Cu,Zn superoxide dismutase, comparing the water molecule sites determined using 1.0 ns MD simulation with those detected by X-ray crystallography. Of the water molecules detected by the two techniques, 20% fall at common sites. These are evenly distributed over the protein surface and located around crevices, which represent the preferred hydration sites. The water mean residence time, estimated by means of a survival probability function on a given protein hydration shell, is relatively short and increases for low accessibility sites constituted by polar atoms. Water molecules trapped in the dimeric protein intersubunit cavity, as identified in the crystal structure, display a trajectory mainly confined within the cavity. The simulation shows that these water molecules are characterized by relatively short residence times, because they continuously change from one site to another within the cavity, thus hinting at the absence of any relationship between spatial and temporal order for solvent molecules in proximity of protein surface.

Active-site copper and zinc ions modulate the quaternary structure of prokaryotic Cu,Zn superoxide dismutase.

The influence of the constitutive metal ions on the equilibrium properties of dimeric Photobacterium leiognathi Cu,Zn superoxide dismutase has been studied for the wild-type and for two mutant protein forms bearing a negative charge in the amino acid clusters at the dimer association interface. Depletion of copper and zinc dissociates the two mutant proteins into monomers, which reassemble toward the dimeric state upon addition of stoichiometric amounts of zinc. Pressure-dependent dissociation is observed for the copper-depleted wild-type and mutated enzymes, as monitored by the fluorescence shift of a unique tryptophan residue located at the subunit association interface. The spectral shift occurs slowly, reaching a plateau after 15-20 minutes, and is fully reversible. The recovery of the original fluorescence properties, after decompression, is fast (less than four minutes), suggesting that the isolated subunit has a relatively stable structure, and excluding the presence of stable intermediates during the dimer-monomer transition. The dimer dissociation process is still incomplete at 6.5 kbar for the copper-depleted wild-type and mutated enzymes, at variance with what is generally observed for oligomeric proteins that dissociate below 3 kbar. Measurement of the degree of dissociation, at two different protein concentrations, allows us to calculate the standard volume variation upon association, Delta V, and the dissociation constant K(d0), at atmospheric pressure, (25 ml/mol and 3 x 10(-7)M, respectively). The holoprotein is fully dimeric even at 6.5 kbar, which allows us to evaluate a lower Delta G degrees limit of 11.5 kcal/mol, corresponding to a dissociation constant K(d0)<10(-9)M.

Local dynamic properties of the heme pocket in native and solvent-induced molten-globule-like states of cytochrome c.

We report the Soret absorption band, down to cryogenic temperature, of native and molten-globule-like state of horse heart cytochrome c. The band profile is analyzed in terms of vibronic coupling of the heme normal modes to the electronic transition in the framework of the Franck-Condon approximation. From the temperature dependence of the Gaussian broadening and of the peak position, we obtain information on the 'bath' of low frequency harmonic motions of the heme group within the heme pocket. The reported data indicate that, compared to the native state, the less rigid tertiary structure of the molten globule is reflected in a higher flexibility of the heme pocket and in greater conformational disorder, allowing the transduction of large-amplitude motion of the protein to the dynamics of the heme pocket.

Flexibility in monomeric Cu,Zn superoxide dismutase detected by limited proteolysis and molecular dynamics simulation.

Limited proteolysis by trypsin of monomeric Cu,Zn superoxide dismutase from Escherichia coli induces a specific cleavage of the polypeptide chain at the level of Lys60 located in the S-S subloop of loop 6,5 where, when compared to the eukaryotic enzyme, a seven-residues insertion, completely exposed to the solvent, is observed. This result suggests that this subloop is disordered and flexible, thus enabling binding and adaptation to the active site of the proteolytic enzyme. Indeed, molecular dynamics simulation indicates that the S-S subloop undergoes high fluctuations and that its high flexibility coupled to an high solvent accessibility can explain the specific bond selection of the protease. As a matter of fact, of the possible 14 solvent accessible proteolytic sites only the Lys60 flexible site is cleaved. High flexibility and solvent exposure are confirmed by the short water residence time for the residues corresponding to the cleavage site evaluated by molecular dynamics simulation. These experiments demonstrate that molecular dynamics simulation and limited proteolysis are complementary and unambiguous tools to identify flexible sites in proteins.

Superefficient enzymes.

Diffusion-controlled enzymes are characterized by second-order rate constants in the range 10(8)-10(10) M(-1)s(-1). These values are at the upper end of the observed rates of many enzyme-substrate reactions and have been predicted by theoretical studies on bimolecular reaction in solution. Such enzymes are considered to be perfect, since their rate-limiting step is not due to any chemical event but to the diffusional association rate between the enzyme and the substrate. Often the enzyme-substrate encounter is facilitated either through the presence of a strong attractive electric field, produced by charges on the enzyme surface, or through the reduction in the dimension of the search process. Here we provide a brief review of some of the enzymes characterized by a very fast second-order constant, focusing attention on triose phosphate isomerase and Cu,Zn superoxide dismutase taken as typical examples of such highly tuned enzymes.